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ABSTRACT 

Our aim is to investigate whether the presence of baryons can have any signif- 
icant influence on the properties of the local Hubble flow which has proved to 
be "cold" . We use two cosmological zoom simulations in the standard ACDM 
cosmology with the same set of initial conditions to study the formation of a 
local group-like system within a sphere of ~ 7 hr 1 Mpc. The first one is a pure 
dark matter simulation (tutidm) while a complete treatment of the physics 
of baryons is introduced in the second one {TUHb aryons ) . A simple algorithm 
based on particles identity allows us to match haloes from the two runs. We 
found that galaxies identified in run& ar2/(ms and their corresponding dark mat- 
ter haloes in tutidm have very similar spatial distributions and dynamical 
properties on large scales. Then, when analyzing the velocity field and the 
deviation from a pure Hubble flow in both simulations, namely when comput- 
ing the dispersion of peculiar velocities of galaxies cr*(i?) and those of their 
corresponding dark matter haloes <jdm{R) in tuudm, wc found no particular 
differences for distances R =1 to 8 Mpc from the local group mass center. This 
suggests that the presence of baryons have no noticeable impact on the global 
dynamical properties of the local Hubble flow within such distances. Then, the 
results indicate that the "true" <r*(i?) values can be estimated from the pure 
dark matter simulation with a mean error of 3 km/s when dark matter haloes 
are selected with maximum circular velocities of V c >30 km/s, corresponding 
to a population of dark matter haloes in run£, aryo „ s that host galaxies. By in- 
vestigating the properties of the Hubble flow at distances R ~0.7 to 3 Mpc, we 
also found that the estimation of the total mass enclosed at the radius of the 
zero- velocity surface Rq, using the spherical infall model adapted to ACDM, 
can be underestimated by at least 50%. 

Key words: Galaxies: Local Group - Galaxies: haloes - Dark matter - 
Methods: N-body simulations 



1 INTRODUCTION 

The study of dynamical and photometric properties of 
galaxies in the local universe represents an ideal frame- 
work to confront predictions of the standard A cold dark 
matter (ACDM) model with observations. In particular, 
an interesting feature of the dynamical properties of the 
local universe is that the Hubble flow is rather "cold", 
e.g. the dispersion in the peculiar velocities within the 
Hubble flow is quite small, namely < 100 km/s (Sandage 
& Tammann 1975; Giraud 1986; Schlegel et al. 1994; 
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Ekholm et al. 2001, Karachentsev et al. 2003; Maccio, 
Governato & Horellou 2005; Tikhonov & Klypin 2009). 

The presence of the dark energy was proposed as 
a possible explanation for the smoothness of the lo- 
cal Hubble flow, first argue by Baryshev et al. (2001), 
Chernin et al. (2001, 2004, 2007) and Teerikorpi et al. 
(2005) and supported by Maccio, Governato & Horellou 
(2005) using a set of N-body simulations. However, Hoff- 
man et al. (2007) compared results from their own sim- 
ulations of CDM and ACDM cosmologies with identical 
parameters, apart from the presence or not of the cos- 
mological constant term. They claimed that no signifi- 
cant differences were noticed in the velocity flow around 
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galaxies having properties similar to those observed in 
the neighborhood of the Milky Way. A similar conclu- 
sion was obtained more recently by Martinez- Vaquero et 
al. (2009) using a range of N-body simulations in differ- 
ent cold dark matter scenarios. They conclude that the 
main dynamical parameter that can affect the coldness 
of the flow is the relative isolation of the Local Group. 
Another approach was proposed by Peirani & de Freitas 
Pacheco (2006, 2008) who derived a velocity-distance re- 
lation by modifying the Lemaitre-Tolman model (TL) 
with the inclusion of a cosmological constant term. They 
found that this new relation, which describes the be- 
havior of the Hubble flow near the central dominant ob- 
jects (<2-3 Mpc), in addition to that derived from the 
"canonical" model (fl m = f ), provide equally acceptable 
fits to the existing available data. As a consequence, no 
robust conclusion about the effects of the cosmologi- 
cal constant on the dynamics of groups could be estab- 
lished. Moreover, Axenides & Perivolaropoulos (2002) 
studied the dark energy effects in the growth of mat- 
ter fluctuations in a flat universe. They concluded that 
the dark energy can indeed cool the local Hubble flow 
but the parameters required for the predicted velocity 
dispersion to match the observed values are out by ob- 
servations that constrain either the present dark energy 
density or the equation of state parameter w(=P x /e x )- 

Other cosmological models were proposed to study 
the local Hubble flow. Dark matter simulations by Gov- 
ernato et al. (f 997) for cosmological models with fi m = 
1 or fi m = 0.3 are, according to these authors, unable 
to produce systems embedded in regions of "cold" flows, 
i.e., with f-D dispersion velocities of approximately 40- 
50 km/s. Recently, Tikhonov et al. (2009a, 2009b) have 
compared the observed spectrum of minivoids in the lo- 
cal volume with the spectrum of minivoids determined 
from the simulations of CDM or WDM models. They 
found that model predictions and observations match 
very well provided that galaxies can only be hosted by 
dark matter haloes with circular velocities greater than 
20 km/s (for WDM) and 35 km/s (for ACDM). They 
have also derived rms deviations from the Hubble flow 
which seem to be consistent with observational values. 

All those past numerical works have used collision- 
less N-body simulations to investigate the puzzle of cold- 
ness of the local Hubble flow and thus, the velocity dis- 
persions have been derived by considering the dynam- 
ical properties of dark matter haloes supposed to host 
galaxies. However, it is nonetheless not obvious that po- 
sitions and velocities of galaxies and their host dark 
matter haloes are identical. A first element of answer 
was given by Weinberg et al. (2008) who have studied 
the subhalo and galaxy populations in a galaxy group 
simulated either by a collisionless cosmological simula- 
tion or a hydrodynamics simulation with the same ini- 
tial conditions. They found that positions and masses of 
large subhaloes are very similar in both runs while they 
can be different for low mass subhaloes whose orbits can 
be more easily modified by the host halo potential. Past 
works also suggest the existence of bias between the star 
and dark matter components (see for instance, Carlberg, 
Couchman & Thomas f990; Zhao, Jing & Borner 2002; 
Sousbie et al. 2008) . In the present work, we study pos- 



sible effects of the presence of baryons in the dynamical 
properties of the local Hubble flow by using high res- 
olution simulations that include most of the relevant 
physical processes that lead to the formation of galax- 
ies inside dark matter haloes. This paper is organized as 
follows: in Sect. 2, we summarize the numerical method- 
ology; in Sect. 3 we present results on the dispersions of 
peculiar velocities around the mean Hubble local flow 
derived from our numerical models and, finally, in Sect. 
4, our main conclusions are given. 

2 NUMERICAL MODELING 

2.1 Initial conditions and simulation 
methodology 

The cosmological simulations analysed in this work use 
the technique of "zoom" (Tormen, Bouchet & White 
1997). We summarize the main steps here. First, we 
have performed a cosmological simulation for a ACDM 
universe using WMAP5 parameters (Komatsu et al. 
2009), namely f2 M = 0.274, tt A = 0.726, tt B = 0.0456, 
H = 70.5 km/s/Mpc, n = 0.96 and a 8 = 0.812. 
The simulation was performed in a periodic box of side 
100 ft -1 Mpc with 512 3 dark matter particles (e.g. with 
mass resolution of ~ 5.6 x 10 s h -1 Mq). The simulation 
started at z ~ 40 and ended at the present time z — 
where we have selected one local group (LG) type halo 
(see section 2.2). In a second step, we have resimulated 
this LG candidate by using an equivalent resolution of 
2048 3 particles within a sphere of 7 ft -1 Mpc. The par- 
ticle mass resolution is progressively degraded around 
this region of interest to reach the lowest resolution of 
256 3 effective particles. Initial conditions have been gen- 
erated from the MPgrafic code (Prunet et al. 2008), a 
parallel (MPI) version of Graf ic (Bertschinger 2001). 

In order to study effects of baryons, we have run 2 
simulations using the same set of initial conditions: in 
the first one, only the dark matter component is con- 
sidered (referred to as runojif) whereas baryons are in- 
cluded in the second one (referred to as runi, ar!/0 „ 3 ) with 
a full treatment of the physics of baryons (see below). 
In these cases, the highest mass resolution are tudm ~ 
8.8xl0 6 ft -I M s (runcaf) &ndm DM « 7.4xl0 6 /i _1 M , 
m aas = m s tars ~ 1-5 x 10 6 h~ 1 M e (run baryon3 ). The 
Plummer-equivalent force softening adopted for the high 
mass resolution particles were 1 and 0.5 kpc for 
dark matter and gas particles respectively and were kept 
constant in comoving units. An intermediate pure dark 
matter zoom simulation with a lower resolution (e.g. 
1024 3 DM effective particles in the central part) has 
been realized in order to check the convergence of the 
dynamical and physical properties of the cental main 
haloes (masses, relative distance of haloes, etc.). 

The simulations were performed using GADGET2 
(Springel 2005). As far as the run(, arl)ons is concerned, 
we have added prescriptions for cooling, star formation, 
feedback from Type la and II supernovae (SN), UV 
background and metal enrichment. Most of this imple- 
mentation have been tested (effects of mass resolution 
on the star formation rate for instance) and used in pre- 
vious works (see, for instance, Peirani et al. 2009a, b). 
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Figure 1. The projected density of dark matter from runuM (left panel) and run(, ar!/ons (middle panel) within a sphere of 
radius 4 Mpc. The right panel represents the projected density of stars in the same volume. 



For the sake of clarity, we summarize the main points 
in the following. 

The evolution of the thermodynamic properties 
of the gas depends on a range of physical processes. 
First we use a standard UV-background prescription 
described in Haardt & Madau (1996). The cooling and 
star formation (SF) recipes follow the prescriptions of 
Thomas & Couchman (1992) and Katz, Weinberg & 
Hernquist (1996), respectively. Gas particles with T > 
10 4 K cool at constant density with a given metallic- 
ity for the duration of a timestep. Gas particles with 
T < 2 x 10 4 A', number density n > 0.1cm -3 , overden- 
sity Apgas > 100, and V.w < form stars according to 
the standard SFR prescription: dp*/dt = c*p ga s/tdyn, 
where p„ refers to the stellar density, tdyn is the dy- 
namical timescale of the gas, and c, is the SF efficiency. 
Instead of creating new (lighter) star particles, we im- 
plement the SF prescription in a probabilistic fashion. 
Assuming a constant, dynamical time scale across the 
timestep, the probability p that a gas particle forms 
stars in a time At is p = 1 — exp(— c» At /tdyn)- For each 
eligible gas particle, we select a random number (r) be- 
tween and 1 and convert it into a star if r < p. 

Next, star particles are treated as evolving stel- 
lar populations that inject energy into the inter-stellar 
medium, through SN explosions, as a function of time. 
For this, we consider stellar lifetimes in the mass ranges 
of 0.8 M e < m < 8.0 M and 8.0 M@ < m < 80.0 M© 
for Type la and Type II progenitors, respectively. Using 
a Salpeter initial mass function for Type II SN, we find 
that the rate of energy injection is: 



H SNla =4.8 x 10- 20 (|^)i5 gJ v(^) 3/2 erg.s-\ (2) 

Equations (1) and (2) are used to compute the en- 
ergy released by SN derived from a star particle i (Ei). 
A fraction 7 of this energy is deposited in the ] th neigh- 
bour gas particle by applying a radial kick to its velocity 
with a magnitude Avj = -J (2wj -yEi /rrij ) , where Wj is 
the weighting based on the smoothing kernel and rrij is 
the mass of gas particle j. We note that all gas neigh- 
bours are located in a sphere of radius Rsn, centered on 
the SN progenitor, to avoid spurious injection of energy 
outside the SN's region of influence. In the following, we 
use the following standard values: 7 = 0.1, Rsn = 0.4 
kpc, and c» = 0.02. 

Supernovae do not only inject energy to the ISM 
but also diffuse heavy metals leading to a progressive en- 
richment of their nearby environment. To estimate the 
stellar metallicity, we assume that the first generation 
of stars will have "zero-metallicity" and the subsequent 
generations will be formed in a medium enriched in met- 
als by the previous generations. To estimate the metal 
abundance, we have considered two elements: Mg (mag- 
nesium), representing the "a-abundance" and Fe (iron), 
representing all heavy elements of the "iron-peak" . The 
former is produced essentially by the type II events 
whereas the later is produced mainly by type la super- 
novae. 

The amounts of Fe and Mg produced by type II su- 
pernovae originating from each stellar cluster are taken 
from Umeda et al. (2002): 



MFe,Mg = (mF e,M 9) g N m*As N , j M®, 



(3) 



HsNj, 



erg.s 



where Esn = 10 erg, m* is the mass of the stellar 
population, and 3.53 < r < 29 Myr. For Type la SN, 
the heating is delayed, since they appear to = 0.8 — 
1.0 Gyr after the onset of star formation. Following de 
Freitas Pacheco (1998), the probability of one event on 
a timescale r after the onset of star formation is given 
by: 



1 ] (1) where m, is the mass of star particle in the simula- 
tion. All Fe (or Mg) is produced after only 29 Myr. 
Using a Salpeter IMF, the fraction by mass of these 
events is \sn ii = 5.72 x 10 _3 Mg 1 and the average Fe 
and Mg masses produced by type II supernovae are < 
m Fe >SNu= 0.068 Mq and < m Mg >sn u = 0.12 M . 
In the case of type la supernovae, the amount of Fe and 
Mg ejected after a given age of a given generation is 



Mfr,Mq — ( m - Fe . A/ 9) gJV m *^SN Ia 1— ( — ) 



'*0\l/2' 



M e , (4) 
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Table 1. Physical properties of the two central main objects in both runs. Numbers in parenthesis correspond to number of 
particles of each component. Dmw—M'H specifies the physical distance between the two central objects and the unit of mass 
is h" 1 M©. 





run DM 




run 6a 


ryons 




Galaxy 


Mtot 


Mtot 


M DM 


M s tar 


a s 


MW 


8.85 x 10 11 


8.96 x 10 11 


7.58 x 10 11 


1.08 x 10 11 


3.02 X 10 10 




(100 097) 


(196 947) 


(102 846) 


(73 600) 


(20 501) 


M31 


8.84 x 10 11 


8.65 x 10 11 


7.43 x 10 11 


8.44 x 10 10 


3.81 x 10 10 




(100 052) 


(183 926) 


(100 856) 


(57223) 


(25 847) 


DMW-M31 


1.01 Mpc 




1.03 


Mpc 





where \sN Ia = 7.2 x 10" 4 Mq 1 , < m Fe >sN Ia = 
0.85 M and < m Mg >SN Ia = 0.026 M© (from model 
W7 by Nomoto et al. 1997). 

The mass of Fe and Mg injected in the medium 
follows the supernovae rate estimated previously, being 
redistributed according to the smoothing kernel among 
the gas particles in the sphere of radius Rsn- Then, 
at each timestep, one can estimate the [Fe/H], [Mg/H] 
and [Mg/Fe] abundances of each gas particle as well 
as a global metallicity Z by using the formula logZ = 
[Fe/H] - 1.70 + log{0.362 + 0.638/ Q ) derived in Salaris, 
Chieffi & Straniero (1993), where f a = 10 [a/Fel is the 
enhancement factor of o-elements with respect to iron. 
In our model, the [Mg/Fe] abundance is used as a proxy 
of [a/Fe]. The inclusion of metal enrichment processes in 
the present study allows us to improve the treatment of 
the gas cooling using metal-dependent cooling functions 
adapted from Surtherland & Dopita (1993). Note that 
a similar methodology was adopted in previous works 
(see for instance, Kobayashi 2004; Tornatore et al. 2004; 
Scannapieco et al. 2005). 

2.2 Properties of simulated Local Group type 
objects 

In this study, we focus on pairs of galaxies with physi- 
cal characteristics similar to the Milky Way- Andromeda 
pair (MW-M31). Each LG group candidate must satisfy 
the following criteria: 

(a) The group contains two main haloes of compara- 
ble mass (~ 0.7—1.2 x 10 12 /i~ 1 Mq), are separated from 
a distance 0.8 < D < 1.0 Mpc and their relative radial 
velocity is negative. 

(b) There is no object with mass greater than 5 x 
W 11 h~ 1 M and 10 13 /i~ 1 M Q within the distance 3 and 
8 Mpc from the center of mass of the main pair of haloes 
respectively. 

(c) There is only one Virgo-like halo (M > 
10 14 l-r 1 M Q ) at the distance of 8 to 13 Mpc. 

Such selection criteria are quite similar to those 
used in previous works (Maccio, Governato & Horellou 
2005; Hoffman et al. 2008; Tikhonov & Klypin 2009; 
Martinez- Vaquero et al. 2007, 2009) and then make the 
comparison of the results easier. Several candidates of 
pairs of dark matter haloes that satisfy criterion (a) can 



be easily found but almost all of them are ruled out by 
criteria (b) and (c). In the end, only one candidate that 
satisfy all of the above criteria has been found. This 
latter system has been resimulated with and without 
baryons and using higher resolutions. The simulations 
were performed on local parallel supercomputers at IAP 
(France) and the CPU times for run^M and runf, arl)ons 
were ~ 2 days and ~ 2.5 months respectively using 
32 processors. Table 1 summarizes the main properties 
of the central objects in both runs and the number of 
particles considered for each component. In the follow- 
ing, we will mostly characterize dark matter haloes by 
their maximum circular velocity which can defined by 
Vmax = max(y / ' Grn(< r)/r), where m(< r) is the mass 
enclosed at radius r. Vmax can be estimated without any 
accurate estimate of the physical boundary of the ob- 
jects which can be proved to be difficult especially for 
subhaloes. 

Fig. 1 shows the projected density of dark matter 
from runoA/ and runbaryons as well as the projected 
density of stars within a sphere of physical radius 4 Mpc. 
The distributions of dark matter look very similar be- 
tween the two simulations. However, the dark matter 
haloes seem to be more concentrated when baryons are 
included in the simulation. This can be explained by 
the fact that dissipation of the gas (from radiative cool- 
ing processes) and subsequent star formation lead to a 
steeper DM density profile due to adiabatic contraction. 
This process can also be observed in dark matter sub- 
structures which become less sensitive to tidal torques 
and will survive a longer time in their host haloes, as has 
been extensively discussed in the past literature (Blu- 
menthal et al. 1986; Gnedin et al. 2004; Duffy et al. 2010 
and references therein). 



2.3 Mock catalogs of dark matter haloes and 
galaxies 

One of the main purposes of this work is to study how 
the inclusion of baryons may affect the dynamical prop- 
erties of the simulated local universe. Then, we have 
used a simple algorithm which allows a direct compari- 
son between the two simulations. 

First, we have prepared mock catalogs of dark mat- 
ter haloes and galaxies at z = for each simulation. To 
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do this we have used the public code AMIGAt (Knoll- 
mann & Knebe 2009) which is an efficient halo and sub- 
structure finder by identifying local density maxima in 
an adaptively smoothed density field. The use of this 
code is also attractive since subhaloes are separated 
form their parent haloes and each (sub)halo center is 
estimated accurately by an iterative algorithm, which 
is crucial in this study. For each object, we define the 
virial radius Rv (and thus the virial mass My) as the 
radius where the enclosed mean density My /(4tt7?^/3) 
is A c times the critical density, p c (z) = 3H(z) 2 /8ttG, 
where H(z) = Ho^jQ, m (l + z) 3 + ^a- A c can be com- 
puted from the spherical top-hat collapse model, and in 
the case of a flat cosmology with f2 m + S1a = 1, one can 
use the fitting formula from Bryan & Norman (1998) 
which depends on both cosmology and redshift: 

A c = 187r 2 + 82a;-39a; 2 , (5) 

x = Q m (z) - 1, (6) 

n m (z) = n m (l + Z ) 3 (-^). (7) 

In the cosmology used for the present study, A c = 97.6 
at z = 0. Then, only bound structures with at least 
100 dark matter particles (e.g. a mass higher than 
8.8 x 10 s /i _1 M Q in runoff) were retained in the dif- 
ferent catalogs. In the following, SPH haloes refers to 
objects identified in ruribaryons which consists of a sin- 
gle dark matter halo containing a fraction of gas and 
(or not) a galaxy, namely a bound structure with at 
least 80 star particles (corresponding to a mass of 
1.2 x 10 s H^Mq). 

Once the two catalogs are built, we match the SPH 
haloes to their equivalent DM haloes of run_DM by us- 
ing the following algorithm. Since each particle in the 
simulation can be identified, it is possible to obtain the 
constitution of each SPH and DM halo at z — 0. If more 
than 50% of dark matter particles of a given SPH halo 
is found in a dark matter halo of tutlom, we assume 
that both objects are the same. To illustrate this, Fig. 2 
shows the distribution of all galaxies that have a corre- 
sponding dark matter halo in runoM within a distance 
of 4 Mpc from the LG mass center. However, it is worth 
mentioning that a small number of SPH haloes have 
no corresponding dark matter halo in tutidm- These 
SPH haloes are mainly substructures belonging to the 
most massive objects such as the MW or M31 haloes. In 
run_DM, they don't exist anymore, mainly because they 
have been completely disrupted by tidal torques. Such 
a statement has already been studied in great detail 
in previous works (see for instance, Nagai & Kravtsov 
2005; Libeskind et al. 2010, Klimentowski et al. 2010). 

It is now interesting to study the statistical differ- 
ences between the radial distances (in the frame of the 
Local Group) and peculiar velocity of galaxies (r g ,Vg) 
identified from runt„ r!(1 , ns and those of their correspond- 
ing dark matter haloes in rmiDM {r dm -Vdm) at z = 0. 
The main results are shown in Fig. 3 where the probabil- 
ity distribution functions of the position ratio tdm /r g , 

t http://popia.ft.uam.es/AMIGA 
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Figure 2. Distributions of galaxies (upper panel) and their 
corresponding dark matter haloes in run dm (lower panel) 
and using the same color code. For clarity reason we have 
limited to a distance of 4 Mpc from the LG mass center. 

the velocity norm ratio Vdm /V g and the angle (V^, Vdm) 
are plotted. Fig. 3 also shows the variations of VDM/r g , 
Vdm/V 9 and (V^,Vda/) with respect to positions r g . It 
appears clearly that most of the galaxies have very simi- 
lar position and velocity properties to their correspond- 
ing dark matter haloes. However, in some cases, we no- 
tice important differences between the two runs. This 
happens mainly for substructures (~ 5 x 1O 9 M0) in 
more massive objects where their dynamics and orbits 
are more likely to be modified by the halo potential, as 
previously observed by Weinberg et al. (2008). More- 
over, a detailed study by Libeskind et al. (2010) has 
also shown that subhaloes in hydrodynamic runs tend 
to be more radially concentrated than those in pure dark 
matter ones. 
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r g [Mpc] Angle (deg) 

Figure 3. Comparison between the radial distances r g (in 
the frame of the Local group) and the peculiar velocities 
V g of galaxies identified in nni(, orvofta and those of their 
corresponding dark matter haloes in run^Af {tdMiVdm) at 
2 = 0: the first column shows the variations of the norm 
ratio roM/fg (first line), the velocity norm ratio Vdm/Vq 
(second line) and angle between velocity vectors (third line) 
as respect to r g while the probability distribution functions 
of tijmI t b ( nrs t line), VnM/Vg (second line) and the angle 
(third line) are plotted in the second column. The red crosses 
correspond to dark matter subhaloes with a mass lower that 
5 x 1O 9 M . 



Among the population of SPH haloes, there is a 
subpopulation of dark matter haloes that contain no 
stars and thus no galaxy. From a theoretical point of 
view, the existence of such dark halo populations can 
be explained by the fact that their star formation have 
been suppressed by feedback processes such as UV- 
background (Thoul & Weinberg 1996; Bullock, Kravtsov 
& Weinberg 2000; Somerville 2002; Ricotti, Gnedin & 
Shull 2002; Benson et al. 2002; Read, Pontzen & Viel 
2006; Hoeft et al. 2006; Okamoto, Gao & Theuns 2008). 
In particular, the presence of a photoionizing back- 
ground suppresses the formation of galaxies with circu- 
lar velocities Vdrc < 30 km/s (Thoul & Weinberg 1996; 
Bullock, Kravtsov & Weinberg 2000) and a similar re- 
sult has been derived by Tikhonov & Klypin (2009) in 
order to explain the observed size of local voids, using 
models with og = 0.75. Moreover, Hoeft et al. (2006) 



have also studied the formation of dwarf galaxies in 
voids using high-resolution hydrodynamical simulations 
including UV-background. They derived a critical mass 
for haloes at z = 0, M c = 6.5 x 10 9 h~ 1 M (s , below which 
the cooling of gas is suppressed and similar results have 
been derived by Okamoto, Gao & Theuns 2008). Then, 
we have compare the properties of dark matter haloes 
of runuM that have either a corresponding galaxies or a 
dark halo in runi, ar!/ons . In Fig. 4, the probability distri- 
bution function is represented for the maximum circular 
velocity V c of DM haloes from tuiidm and the variations 
of V c as respect to their total halo mass. The results sug- 
gest that there is a critical value [V c ~ 30 km/s) where 
the star formation is suppressed, in good agreement with 
past works. Moreover, dark matter haloes that have cor- 
responding starless SPH haloes 111 VMUbaryons <W6 char- 
acterized by Vc < 35 km/s and have a mass lower than 
~ 7 x 10 9 h~ x Mq. Note that, the circular velocity of 
dark matter haloes of mnoM can be well approximated 
by V c = 39.08 km s" 1 x (M/10 10 h- 1 M @ ) - 32 in good 
agreement with Hoeft et al. (2006) and Klypin, Trujillo- 
Gomez & Primack (2010). 



3 THE LOCAL HUBBLE FLOW 

The local Hubble flow has proved to be cold according 
to past observational analysis. As mentioned in the in- 
troduction, this puzzle has been studied in various nu- 
merical works but no one so far has investigated the 
impact of baryons, the main purpose of this work. 

In the following, the peculiar velocity dispersion oh 
in the local volume is estimated by using the root mean 
square (RMS) of the residual from the global Hubble 
flow, defined by: 

1 N 

oh = ^2 (vi - v H (ri)) (8) 

i=l 

where Vi and n are radial velocities and radial dis- 
tances of galaxies (or dark matter haloes) in the frame 
of the Local Group where the origin matches with the 
mass center of the MW-M31 pair, and vh is the global 
Hubble flow. In order to compare with recent works, 
and more specifically with the analysis of Tikhonov & 
Klypin (2009) and Martinez- Vaquero et al. (2009), we 
have estimated oh using different samples of dark mat- 
ter haloes or galaxies with distances between 1 to D Mpc 
(where D=2, 3, 4, 5, 6, 7 and 8 Mpc) and considering 
vn(ri) = HoVi, with Ho = 70.5 km/s/Mpc. However, it 
is worth mentioning that for distances within ~ 3 Mpc, 
due to the decelerating influence of the main central ob- 
jects to nearby low mass objects, the Hubble flow is not 
well described by vn(ri) = Hort (see Fig. 7) and one 
should use another relation such as the one derived in 
Peirani & de Freitas Pacheco (2008). Then, the study 
has been limited to a distance of 8 Mpc to avoid to be 
in the vicinity of the boundary of the highest resolution 
domain. First, sample A includes all galaxies identified 
in runt aryons and represents therefore the fiducial sam- 
ple. For this case only, oh has been estimated either 
by assuming a 10% error as the mean error of distance 
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Table 2. The velocity dispersion around the Hubble flow vn(r) = Hor for different samples of selected galaxies (runj,„ rBO „ s ) or 
dark matter haloes (runrjM) as a function of the distance from the mass center of the simulated MW-M31 pairs. Each number 
in parenthesis gives the number of objects that have been considered in the computation of an ■ 
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a cr H computed from galaxies identified in nm(, OT . !/tm . a (Sample A). 

b &[f computed from galaxies identified in runi, aryons (Sample A) but corrected from distance errors. 
°<IB computed from galaxies identified in runt, ar j /ons that have a corresponding DM halo in runrjM (Sample B\). 
d au computed from DM haloes identified in runrjM that have a corresponding galaxy in run f)ar!/onlJ (Sample B^)- 
e tJH computed from DM haloes identified in runrjM with a maximum circular velocity V c >45 km/s (Sample Ci). 
' uh computed from DM haloes identified in rxirioM with a maximum circular velocity Vc >35 km/s (Sample C-i). 
B rJH computed from DM haloes identified in runrjM with a maximum circular velocity V c >20 km/s (Sample C3). 

h rms = ^EKM- ff 'M) 



measurements (referred to as cr*, cor ) or with no error 
(referred to as o,). Samples B\ and B2 regroup galaxies 
and their corresponding dark matter haloes in rtuiDM 
respectively and have therefore the same number of ob- 
jects which permits a direct comparison between the 
two runs. Finally, samples C\, C2 and C3 gather dark 
matter haloes extracted from rxmoM with circular ve- 
locities higher than 45, 30 and 20 km/s respectively. 

Fig. 5 shows the Hubble diagram for galaxies and 
dark matter haloes belonging to sample A, Bi and £>2. 
As expected from results presented in the previous sec- 
tion, radial distances and velocities of galaxies and those 
of their corresponding dark matter haloes in runoM fol- 
low the same trend. Such a statement is also confirmed 
by Table 2 where all the oh values derived at different 
distances from the LG mass center and for different sam- 
ples are shown. In particular, oh values derived from 
samples Bi and B2 are very close and suggests that the 
presence of baryons does not affect the velocities dis- 
persion of satellites around the main pair at distances 
1-8 Mpc. Moreover, results from samples Ci, C2 and C3 
show that estimations of oh from the pure dark matter 
simulation give similar values than those derived from 
the fiducial sample as it can been clearly seen in Fig. 6. 
For instance, when dark haloes are selected with V c > 30 
km/s, the number of objects considered is quite close to 
the "true" number of galaxies and the mean error on 
oh is in order of ~ 3 km/s. However, when haloes are 
selected with V c > 45 km/s and V c > 20 km/s, both 
cases lead to an estimation of oh with a mean error of 
~ 5 km/s as respect to <r* values although the number 
of objects taken into account for the oh computation 
are underestimated or overestimated respectively. 

To finish, it is instructive to compare the veloc- 
ity dispersion values derived from our numerical model 
to oh values obtained from the analysis by Tikhonov 
& Klypin (2009) who have considered the most com- 
plete observational data using an updated version of 



the Karachentsev et al. (2004) sample of galaxies. The 
number of galaxies used to compute o, (first column 
of Table 2) is in quite good agreement with number of 
galaxies considered within 1-8 Mpc (see their Table. 3). 
However, we found that oh values derived in this work 
tend to be (~ 10%) smaller within 1-8 Mpc. We have 
also compared in Fig. 7 our simulated data to updated 
data on distances and velocities of galaxies for distances 
0.7 < D < 3 Mpc without an additional dipole compo- 
nent (Karachentsev et al. 2009). The nice agreement 
between simulated data and observations lead to simi- 
lar values of the dispersion of peculiar velocities. Indeed, 
using in this case vn(r) = V r (r) derived in Peirani & 
de Freitas Pacheco (2008), we found oh = 29.3 km/s 
from observational data whereas a* = 33.0 km/s and 
o"dm,v c >30 = 29.9 km/s respectively for galaxies from 
and dark matter haloes from runrjM with 
V c > 30 km/s. But an important feature of the Hubble 
flow at such distances is the decelerating influence of the 
central main objects on nearby galaxies which lead the 
radial velocities-distances regression to cross the line of 
zero velocity at the turnaround radius Ro. In this re- 
gard, Lynden-Bell (1981) and Sandage (1986) proposed 
an alternative method to the virial relation in order to 
estimate the mass of the Local Group, which can be 
extended to other systems dominated either by one or 
a pair of galaxies. Based on the Lemaitre-Tolman (LT) 
model (Lemaitre 1933; Tolman 1934), they show that, if 
the velocity field close to the main central body (probed 
by satellites) allows the determination of Ro, then the 
total central mass can be calculated straightforwardly 
from the relation 



where To is the age of the universe and G is the grav- 
itational constant. In a flat cosmological model with a 
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Figure 4. Dark matter haloes from runr)M : the probability 
distribution functions of their maximum circular velocity V c 
(upper panel) and the variations of their mass as respect to 
V c (lower panel). Blue color refers to DM haloes that have 
corresponding starless DM haloes in runi, ar . yo „ s whereas red 
color refers to DM haloes that match SPH haloes with a 
galaxy. In the lower panel, the solid line represents our best 
fit (see the text for more details) while the dashed lines show 
V c = 35 km/s and M = 7 X 10 9 h.- 1 M Q . 



cosmological constant A, equation (9) needs to be mod- 
ified as follows (Karachentsev et al. 2007): 



M 



TT 2 RqHq 

8GT 2 /(Oy 



where 

/(12a) = o ^rrrcosh 

12a 



2<> 



3/2 



- i 



(10) 



(ii) 



We have estimated Ro from our two mock catalogs 
of galaxies and haloes by fitting the V r = V r (r) formula 
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Figure 5. Simulated radial velocity-distance data for galax- 
ies of sample B\ (green dots) and dark matter haloes from 
sample B2 (red dots). The blue dots correspond to velocity- 
distance data for galaxies that don't have a corresponding 
halo in runrjM- The solid line represents the Hot Hubble 
flow and the dashed ones show -Hor ± 98.9 km/s. 



derived in Peirani & de Freitas Pacheco (2008) to our 
simulated data. From the best fit solutions, we obtained 
Ro = 0.98 Mpc (galaxies sample) and Rq = 1.03 Mpc 
(for dark matter haloes sample with V c > 30 km/s). 
Using Qm = 0.274 and Ho = 70.5 km/s/Mpc, one can 
derive from Eq. (10) the total masses enclosed at radius 
Ro, namely M = 2.0 x 10 12 M Q (galaxy sample) and 
M = 2.3 x 10 12 Mq (dark matter haloes sample) respec- 
tively. However, the total mass enclosed in a radius of 
1 Mpc is ~ 3.4 x 10 12 M Q from both simulations and 
suggest that masses derived from Eq. (10) are underes- 
timated by a factor 1.5-1.7. 



4 DISCUSSION AND CONCLUSIONS 

In the present paper, we have used cosmological N-body 
simulations with and without a complete treatment of 
the physics of baryons to study the formation of a local 
group type halo. We have first identified a LG system in 
a simulation with a low mass particle resolution, using 
standard criteria. Then, the best candidate has been 
re-simulated with higher resolution using the technique 
of zoom. This approach differs from the technique of 
constrained simulations used for instance in the CLUES 
project (Yepes et al. 2009), which include the correct 
motion and position of objects in a large volume such as 
the Local Supercluster, the Virgo and the Coma cluster 
and the Great Attractor (see also Lavaux 2010). 

We focus on the dynamical properties of the local 
Hubble flow and in particular, the influence of baryons. 
From the two runs, we built mock catalogs of either 
galaxies or dark matter haloes and computed their ve- 
locities as respect to the mass center of the central main 
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Figure 6. The velocity dispersion around the Hubble flow 
vh(t) = Hor as a function of the distance from the mass 
center of simulated MW-M31 pairs for galaxies identified in 
YVLnt, aryona (o*, black line), and for dark matter haloes from 
runrjAf with a maximum circular velocity defined by V c > 20 
km/s (blue circles), V c > 30 km/s (green circles), V c > 45 
km/s (red circles). The dashed lines show <t„ ± 5 km/s. 



pair (Milky Way-M31). We found that the dispersions of 
peculiar velocities around the mean local Hubble flow of 
galaxies in run6„ r!/0 „ s and those of their corresponding- 
dark matter haloes in run_oj\/ are very close for dis- 
tances D = 1 to 8 Mpc. Moreover, similar <jh values 
(with mean errors < 5 km/s) are also obtained for sam- 
ples of dark matter haloes extracted from runrjA/ with 
constraint of their maximum circular velocity. These re- 
sults suggest that the global dynamical properties of the 
Hubble flow is not affected by the presence of baryons. 
Such a result is not very surprising since galaxies are 
expected to form in the gravitational potential well of 
dark matter haloes. Therefore, their properties are ex- 
pected to be closely related to those of their host haloes. 
This statement is confirmed on large scales as shown in 
Figs. 2 and 3. However, some discrepancies between po- 
sitions and dynamical properties of galaxies and those of 
their corresponding haloes may raised in the case of low 
mass objects in denser environment. For instance, the 
presence of baryons is supposed to favor the survival of 
substructures inside massive haloes against tidal strip- 
ing (see Duffy et al. 2010 and references therein) and can 
affect their distribution (Weinberg et al. 2008; Libeskind 
et al. 2010, Knebe et al. 2010). But these effects seem 
to have no particular consequences in the estimation of 
the global oh values. It should also be pointed out that 
runi, ar!/ona does not include any prescriptions of strong 
feedback such as galactic winds or Sedov shock waves 
that are expected to lead to a significant ejection of gas 
from the central region of forming galaxies. As a result, 
the fraction of baryons (in the form of stars and cold 
gas) in the central region of simulated massive galaxies 
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Figure 7. Upper panel: radial velocities and distances for 
observed (green circles) and simulated (red circles) galaxies 
in the neighborhood of the LG with respect to its centroid. In 
the lower panel, the velocity and distance data derived from 
dark matter haloes in runjjM with V c > 30 km/s is shown. 
In both panels, the inclined dashed line represents the linear 
Hubble relation with Ho = 70.5 km/s while the solid line 
represents the best fit to the v r = v r (r) relation derived in 
Peirani & de Freitas Pacheco (2008). 



such as the Milky Way or M31 is close to the cosmologi- 
cal baryonic fraction (see Table 1) whereas observations 
suggest lowest values (Fukugita, Hogan & Peebles 1998; 
Bell et al. 2003). However, no one has solved this miss- 
ing baryon problem yet and recent studies even suggest 
that theoretical sources of strong winds such as AGN 
and supernovae are insufficient to explain the missing 
baryons in Milky Way type galaxies (see for instance 
Anderson & Bregman 2010; Silk & Nusser 2010). Al- 
though the present study mainly focus on lower mass 
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galaxies (< lO lo Af0) in which supernovae feedback is 
supposed to have a significant impact in their formation, 
the other sources of strong feedback should be taken 
into account as well. Then, a natural extension of this 
work would be to add these processes while increasing 
the mass resolution in the simulation (see for instance 
Ceverino & Klypin 2009) to see whether the conclusions 
presented here (especially the distribution of galaxies as 
respect to their corresponding dark matter haloes) are 
significantly affected of not. 

In good agreement with previous theoretical and 
numerical works, we found that the halo population in 
the pure dark matter run can be divided into two groups 
according to their maximum circular velocity V c . There 
is indeed a critical value, V c ~ 30—35 km/s (correspond- 
ing to haloes with mass lower than ~ 7 x 10 9 /i _1 Mq 
in run_DM), below which the corresponding SPH haloes 
don't host any galaxy. Such a phenomenon may be ex- 
plained by the fact that the cooling of gas is expected 
to be suppressed by feedback processes such as a UV- 
photoionisation in low mass haloes (see references in 
paragraph 2.3). As regards this point, we cannot ex- 
clude the fact that the numerical resolution used in this 
work may affect the results especially for low mass sys- 
tems. However, haloes with mass of - 7 x 10 9 /i _1 M 
are composed by ~ 800 particles and this seems reason- 
able for the estimation of the thermodynamic properties 
of the gas component. Then, table 2 indicates that the 
number of dark matter haloes in rutin m selected with 
V c > 30 km/s is actually quite close to the real number 
of galaxies selected Tunt, aryons . Therefore, one can use 
this single criteria to select dark matter haloes in a col- 
lisionless simulation to have an good estimation of an- 
The existence of a dark halo population in the Local 
Group can also have important consequences in under- 
standing the discrepancies between the large number 
of subhaloes present in simulations but not observed 
(Kauffmann, White & Guiderdoni 1993; Moore et al. 
1999; Klypin et al. 1999). This problem is beyond the 
scope of this paper and will be investigated in detail in 
a forthcoming paper. 

The velocity dispersion obtained from our fidu- 
cial model (e.g. from galaxies) is in quite good agree- 
ment with observational expectations. For instance, 
with added corrections for distance errors, we obtained 
o",, CO r = 90.4 for D = 1 — 8 Mpc, a value which tend 
to be even lower that the one derived from observa- 
tional analysis, namely 99.3 km/s (Tikhonov & Klypin 
2009) and with a reasonable number of galaxies consid- 
ered. For distances from 0.7 to 3 Mpc, we also found 
a nice agreement between our simulated data and ob- 
servational ones from Karachentsev et al. (2009). Then, 
these results suggest that there is no particular prob- 
lem with the ACDM model in reproducing the veloc- 
ity dispersion derived from observational analysis. A 
similar conclusion was obtained from Martinez- Vaquero 
et al. (2009) who found that a non-negligible fraction 
of LG like objects simulated form various cosmological 
models present a an value close to (and even smaller) 
than the observed value. Moreover, similar results have 
been also obtained from the analysis of Tikhonov & 
Klypin (2009). Since the present study has been lim- 



ited to only one realization, due to the high computa- 
tional cost when the physics of baryons is included in the 
simulations, general conclusions cannot be drawn. How- 
ever, one of our main criteria to select LG candidates 
was their relative isolation (no massive galaxies within 
3 Mpc). We stress that we found only one candidate 
(among 16 MW-M31 pair candidates) that was satis- 
fying all the criteria. This raised the question whether 
the local group is located in a particular place of the 
universe or not. 

To finish, we have also compared our simulated data 
to observational ones for R — 0.7 — 3 Mpc. The nice 
agreement between them allows us to test models, in 
particular the estimation of the total mass enclosed at 
the zero- velocity surface radius Ro using the Lemaitre- 
Tolman model adapted to ACDM. It appears that the 
estimation of the LG mass using this approach is under- 
estimated by at least 50%. This is probably due to the 
fact that, on the one hand, the hypothesis of a spher- 
ical infall collapse is not valid anymore in this specific 
case. On the other hand, the mass enclosed within Rq is 
not constant in the time, as assumed by the TL model 
(e.g. no mass accretion), and this effect may have an 
important influence in the estimation of the final total 
mass. 
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